Conductance oscillations of a spin-orbit stripe with polarized contacts 
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We investigate the linear conductance of a stripe of spin-orbit interaction in a 2D electron gas; 
that is, a 2D region of length I along the transport direction and infinite in the transverse one in 
which a spin-orbit interaction of Rashba type is present. Polarization in the contacts is described 
by means of Zeeman fields. Our model predicts two types of conductance oscillations: Ramsauer 
oscillations in the minority spin transmission, when both spins can propagate, and Fano oscillations 
when only one spin propagates. The latter are due to the spin-orbit coupling with quasibound 
states of the non propagating spin. In the case of polarized contacts in antiparallel configuration 
Fano-like oscillations of the conductance are still made possible by the spin orbit coupling, even 
though no spin component is bound by the contacts. To describe these behaviors we propose a 
simplified model based on an ansatz wave function. In general, we find that the contribution for 
vanishing transverse momentum dominates and defines the conductance oscillations. Regarding the 
oscillations with Rashba coupling intensity, our model confirms the spin transistor behavior, but 
only for high degrees of polarization. Including a position dependent effective mass yields additional 
oscillations due to the mass jumps at the interfaces. 

PACS numbers: 71.70.Ej, 72.25.Dc, 73.63.Nm 



I. INTRODUCTION 

Spintronics attempts to manipulate the electron spin, 
alone or in combination with the electron charge, to tune 
the current in a device or as the bit of information.^ 
This novel technology may lead to faster responses and 
lower power consumptions as compared to the more con- 
ventional electronics. A promising approach to spin- 
tronics is the control of the spin-orbit interaction, par- 
ticularly, of the Rashba interactioni^ This is a spin- 
orbit couphng due to the lack of inversion symmetry in 
semiconductor heterostructures, such as those based on 
InAs or GaAs semiconductors, for which tunability using 
gate contacts has been demonstrated.'^''* Exploiting this 
tunability, Datta and Das proposed the spin field-effect 
transistor)^ a ballistic channel between two ferromagnetic 
leads where current can be manipulated by changing the 
Rashba strength via an external gate. Despite its sim- 
plicity, the limitations of the physical system^ and of the 
ballistic transport- proved to be long-time obstacles to 
realizing this device. A recent experiment, however, has 
overcome these problems demonstrating the feasibility of 
the Datta-Das transistor.-^ 

The Rashba Hamiltonian is composed of two spin- 
dependent terms; one conserves the longitudinal momen- 
tum while the other couples longitudinal and transverse 
momenta. They are usually called precession and mixing 
terms, respectively. The mixing term is not present in 
ideal ID Datta-Das transistors. More realistic quasi-lD 
systems are usually considered including a confinement 
in the transverse direction to transport. The quantum 
wire with homogeneous spin orbit interaction has been 
discussed, for instance, in Refs. I9l4ll[ In a similar quan- 
tum wire configuration but with the Rashba coupling re- 
stricted to a finite region of the wire, recent works showed 



the importance of the mixing term in the modulation 
of the conductance.— "-i^ These modulations are exam- 
ples of Fano resonances due to the coupling with quasi 
bound states. Since this behavior is caused by the spin- 
orbit coupling alone it has been called the Fano-Rashba 
eSectM 

The lateral dimension of the transport channel in the 
experiments by Koo et al^ was several microns, which 
indicates a high degree of 2D character, thus deviating 
from the ID or quasi-lD regimes. The quasi-lD multi- 
channel case with polarized contacts was considered in 
Refs. [iTHigl while the 2D system, without confinement in 
the transverse direction, has been also addressed in Refs. 
l20l - [25l In most cases sharp transitions between the con- 
tacts and the channel are assumed and matching of the 
wave functions at the interfaces is the required condition. 
As an alternative, our approach assumes smooth transi- 
tions and describes transmission and reflection between 
contacts and channel from the numerical wave function, 
solution of the complete Schrodinger equation. Our pur- 
pose is to provide additional insight on the origin and 
characteristics of different types of conductance oscilla- 
tions. 

In this work we extend the analysis of Ref. [191 to the 
case of vanishing transverse confinement. We thus focus 
our interest in a two dimensional electron gas (2DEG) 
with a stripe of spin-orbit interaction. Polarization in 
the contacts is modeled by means of effective Zeeman 
fields, treating the cases of parallel and antiparallel po- 
larizations in arbitrary directions. We assume these Zee- 
man fields in the contacts are the stray fields of nearby 
ferromagnets deposited on top of the 2DEGi-? The di- 
rection of the Zeeman field is modeling the orientation of 
the ferromagnet for each contact. In particular, we are 
interested in parallel and antiparallel orientation of the 



ferromagnets in the plane of the nanostructure. The hn- 
ear conductance of our model system displays two kinds 
of oscillations: Ramsauer oscillations when the two spins 
are propagating and Fano oscillations when only one spin 
can propagate while the other one is evanescent. Look- 
ing at the evolution of the linear conductance with the 
Rashba coupling intensity, our model confirms the Datta- 
Das oscillation of the conductance, but only for high de- 
grees of polarization in the contacts. 

This paper is organized as follows. In Sec. II we de- 
scribe the physical system and the model for the current 
and linear conductance. Section III discusses the depen- 
dence of the conductance on energy while the Datta-Das 
transistor configuration is studied in Sec. IV. In Sec. V 
we study how the results are affected by a position de- 
pendent effective mass. Finally, Sec. VI presents our con- 
clusions. 



II. PHYSICAL SYSTEM AND MODEL 

We consider a semiconductor 2DEG in the xy-plane 
with a region of Rashba spin-orbit interaction shaped like 
an infinite stripe of width £ oriented along y. Figure [T^ 
shows a sketch of the physical system. Transport is along 
X and the asymptotic leads (contacts) are assumed to be 
spin polarized along a given direction h. The system 
Hamiltonian reads 



n = - 



2mQ 

+ A(a;) n ■ a 



Aix)\ +nR, 



where "Hr is the Rashba Hamiltonian, 



"Hfl = - ia{x)pyax - -{a{x),px}ay 



(1) 



(2) 



Polarized leads in the direction of h are described by 
means of a Zeeman field A (a;) that couples to the spin 
vector a. A positive scalar potential |A(2:)| is also intro- 
duced in order to align the majority spin potentials in the 
contacts with the potential bottom of the central region. 
This eliminates the effect of a potential mismatch for this 
spin component and, in practice, it would correspond to 
use a potential gating of the central region. In Eq. ^ 
the functions determining the Hamiltonian are A(a;) and 
a{x). These quantities take a constant value in the three 
parts of our system: left contact (L), central region and 
right contact (R), and they vary smoothly, described by 
a Fermi-type function, at the interfaces. See Appendix 
|X]for the precise definitions. 

We denote by mo the conduction-band effective mass 
of the semiconductor and by ao the Rashba intensity 
of the central region. The Zeeman field in contact c, 
where c = L,R, is denoted by Ac, respectively. The 
case of parallel polarized contacts (P) corresponds to 
Ai = A/{ = Ag, while the case of antiparallel polar- 
izations (AP) corresponds to A^ = — Aj^ = Aq, where 
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FIG. 1: (Color online) Sketch of the physical system (a) and 
of the spatial variation of Rashba intensity a(x) (b) and of 
the spin-dependent potentials ^±(2;) (c,d). See Sec. IL 



Aq is half of the absolute Zeeman splitting. For simplic- 
ity, Aq is assumed equal in both contacts. We use the 
notation nP and nAP to indicate parallel and antiparal- 
lel configurations along a certain direction fi. Figure [T|d 
shows the variation of the Rashba intensity a{x). It also 
shows the potentials Vs, for s = ± spins, defined as 



v,{x) =sA(a;)-t-|A(x)| 



(3) 



Notice that in the P configuration the s — — spin sees 
no potential at all while s = -I- is confined by a potential 
well of width d. On the contrary, in the AP configuration 
both spins feel a potential step, but in opposite contacts. 
As we will discuss below, these differences in potential 
landscape for -I- and — spins greatly infiuence the trans- 
port properties of the stripe with polarized contacts. 

Typical values of the spin-orbit intensity for InAs based 
semiconductors can be tuned around ao ~ lOmeVnm, 
with about one order of magnitude range. Assuming a 
system length of ^ « 1 /xm and a Zeeman splitting of 
Aq ~ 0.3 meV this implies that, in adimensional units, 

one has ao ~ 0.3-W Ti^Ao/rrio and ^ sa 10 yfi^ /nioAo. 
Having in mind these orders of magnitude we shall ex- 
plore the variation with energy, for energies around the 
Zeeman gap, and for Rashba coupling intensities around 
0.5 in adimensional units. 

Asymptotically, in the 2D contacts, the Hamiltonian 
eigenfunctions factorize as a plane wave times a spinor in 
the direction of h, 



<^cs{r,'q;Kcs) = exp (i^cs • r)Xs(??) , 



(4) 



where c = L,R and s = ± are labelling contact and 
spin, respectively. The wavenumber Kcs = {kcs,<lcs) is 
composed of the longitudinal (fc^s) and transverse (qcs) 
components. Anticipating a result emphasized below, we 
note that the transverse momentum is a good quantum 



number of the system Hamiltonian, Eq. ([T]). Therefore, 
Qcs must be a characteristic of the wave function not only 
in an asymptotic region c and for a given spin s, but 
throughout the system, i.e., qcs = q- At a given Fermi 
energy E we then have 



= fc 



q' = 



2mo 



(^-sAc-|Ae|). 



(5) 



The physically acceptable wave functions fulfill 
Schrodinger's equation 



{H-E)^ = 0. 



(6) 



A most general wave function can be taken as a sum on 
spins and over all transverse momenta (an integral in q) 

*(^ '^) = Y. j dq^,s{x) e^'^y xsiv) , (7) 



s=± 



where the unknown functions ipqs{x) can be interpreted 
as the wave amplitudes in each channel given by {qs). 
Projecting Eq. ([5]) we obtain the channel amplitude equa- 
tions 

+ X! 1 ( "(^)'?(ska;|s') + -a'{x){s\ay\s') j ipqs'ix) 



s'=± 



dx 



+ ia{x){s\ay\s')—tPqs'{x) \ = 0. (8) 



Notice that the channel equations for different q's are 
uncoupled due to the translational invariance of the sys- 
tem in the transverse direction. At a given q, however, 
the two spin components do couple with each other due 
to the Rashba spin-orbit interaction. This coupling is 
described in Eq. ([8|) by the matrix elements {s\cFx,y\s') 
which can not be diagonalized simultaneously. For any 
orientation of the spin quantization axis h, therefore, 
there is a Rashba-induced interference of the two spin 
projections. In the contacts the Rashba coupling vanishes 
and the wave function recovers the good spin eigenstates 
<^cs given in Eq. (g]). 

Integration of Eq. (|5]) determines the transmission T^/s, 
which represents an electron entering the system from the 
left contact with spin s and going to the right lead with 
spin s' . It also gives T^',^, i.e., from the right contact 
with spin s to the left one with spin s' . In terms of these 
transmissions the total current Ix, per unit of length in 
the transverse direction Ly, can be obtained by adding 
up the contributions of all electrons in each contact,— 



Ix _ e 
Ly - (2^)2 



jiyi/^^^Y.fl.\^^)n'M. (9) 



In Eq. (HJ , /is (k) represents the distribution function 
of electrons in contact c with spin s, with the upper index 
indicating right (-I-) or left (— ) direction of motion of the 
corresponding electron. The distribution of electrons in 
each contact is given by a Fermi function, characterized 
by a given chemical potential iic- In the linear response 
regime the bias SV = fin — fiL is very small and it is 
enough to retain the linear conductance Ix = GSV . From 
Eq. Q we find the conductance per unit of transverse 
length 



Gr GrQ 

L„ ~ An 



V <^ KLs / d9\ COS eiT.^siuLs:! 



/.3V 2 

+ KRs / de\cose\T;^,^{KRs,( 



(10) 



where Gq = e^ /h is the conductance quantum and k^s is 
the Fermi wavevector in contact c, given by Eq. ([S]) when 
E is the corresponding Fermi energy. 



III. RESULTS 
A. Numerical energy dependence 

Figures [2] and [3] show typical results obtained numer- 
ically from Eq. (|10|) for x- and y-polarized contacts, re- 
spectively. The case of polarization along z is similar to 
the X one. For each direction {x and y) we also compare 
the situation of parallel and antiparallel Zeeman fields in 
the contacts. 

For a vanishing Rashba field the transmissions can be 
obtained analytically. If, in addition to uq = 0, the Zee- 
man fields also vanish (Ao — 0), the transmissions triv- 
ially become one and the exact conductance is then 



^ ^ ^^2 l2m^E 



(11) 



When Aq ^ we have to distinguish P and AP configu- 
rations. The P case is characterized by a perfect trans- 
mission of the — spin, while the -I- spin feels the f + po- 
tential of Fig. [TJ Therefore, its transmission switches on 
only when E > 2Ao. When this occurs, the underly- 
ing potential well makes the transmission of the -I- spin 
oscillate with energy, even with vanishing spin-orbit. Fol- 
lowing Ref. [13, we call these variations Ramsauer oscilla- 
tions, in analogy with the Ramsauer effect in electron 
scattering. 2* The importance of these oscillations was 
pointed out in Ref. [20. Notice also that with vanish- 
ing Rashba field the results for x and y orientations of 
the Zeeman fields are identical. 

The energy 2Ao signals the transition threshold from 
only one propagating spin when E < 2Ao, to both spins 
when E > 2Ao. At ao = Oj the P conductance below 
threshold is given by a pure square root behavior, as 
in Eq. PTj) . while above threshold it shows Ramsauer 
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FIG. 2: (Color online) Conductance as a function of energy for 
polarized contacts along a: in P (upper) and AP (lower) config- 
urations. Gray (red in color) and black curves are the results 
with and without Rashba stripe, respectively. The Zeeman 
field parameter Ao is taken as energy unit and, accordingly, 
Lu = \/?i^/moAo as length unit. The unit of Rashba cou- 
pling intensity is given by au ~ \/h'^Ao/mo. Other Hamil- 
tonian parameters: d — 20 Lu, i ~ 8Lu, a ~ 0.1 Lu- 



oscillations of the minority spin transmission. The AP 
transmission is exactly zero below threshold (ao — 0) 
and above it begins to increase smoothly, as expected for 
a spin valve. Note that the Ramsauer effect is not active 
in the AP configuration since the underlying potential is 
a step, instead of a well. 

Turning now to the spin orbit effects, the most conspic- 
uous one is that for E < 2Ao, when the contacts are fully 
polarized, the Rashba field induces the appearance of os- 
cillations in the xP and xAP configurations (Fig. [2]). As 
we will discuss in detail in the next subsection, these os- 
cillations are due to resonant Fano interferences between 
the propagating spin and the quasibound states of the 
opposite evanescent spin. They are qualitatively similar 
to the Fano- Rashba interferences discussed in Ref. [ij for 
quantum wires. Here, however, the quasibound states are 
caused by the polarized contacts and not by the Rashba 
field itself. In the y orientation (Fig. [S]) , the Fano oscil- 
lations below threshold are absent, and only some small 
variation from the vanishing spin-orbit case can be seen. 
In general, as shown in Figs. [5] and [31 for the P configu- 
ration the results with Rashba coupling (gray-red curve) 
are slightly below the results without spin-orbit (black 
curve); while For the AP configuration the situation is 
reversed. 



FIG. 3: (Color online) Same as Fig. [2] for polarization of the 
contacts along y. 



B. The ansatz model 

Let us focus our attention on the oscillations that ap- 
pear in the fully-polarized-current case, when E < 2Ao, 
both in parallel and antiparallel configurations. The up- 
per panel of Fig. 2] shows the transmission for vanishing 
transverse momentum in xP configuration, in compari- 
son with the total transmission obtained by integrating 
over angle 9 in Eq. pUj) . We clearly see that the oscil- 
lations in G are due to the deep minima in the trans- 
mission for q = 0. Besides, the position of these minima 
coincide with the energies of the bound states in the v+ 
potential (dashed lines). At finite q's, not shown in the 
figure, the transmission minima are shifted or they can 
even disappear. Physically, we indeed expect the q = 
contribution to dominate the conductance since in this 
case all the available energy is used in the longitudinal 
wavenumber. The present transmission minima are ex- 
amples of Fano resonances due to the interference with 
quasibound states. To better understand this behavior 
this subsection presents a simplified model involving the 
quasibound states in an explicit way. 

Assuming q ~ and xP configuration Eq. ([5]) trans- 
forms to 



2?7io dx'^ 



-E\i>^ = V„rij+ , 



(12) 



2mo dx^ 



E] ^, 



-K>^- , (13) 



where we have defined the gradient-dependent mixing po- 
tential V,n = a'{x)/2 + a{x)d/dx. Equations p^ and 
(1131) constitute a two-channel model, where ^_ is prop- 
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FIG. 4: (Color online) Contribution with q = to the 
linear conductance G/Ly. The vertical dashed lines signal 
the bound-state energies of the v+ potential. Except for 
ao = 0.1 au, we have used the same parameters and units 
of Fig. [2] Upper panel: 5 = conductance from Eq. (I10|l . 
For comparison, the gray (red) curve shows the full G/Ly of 
Eq. (|10|l . Lower panel: q — conductance from the ansatz 
model. Solid and dashed lines are from Eqs. (|2U|I and (|21|l . 
respectively. 



agating while "0+ is evanescent, with a localized mixing 
described by Vm- Similar models were obtained for im- 
purities in quantum wires, where semianalytical solutions 
were worked out using Green functions^^ or the ansatz 
by Nockel and Stone.— 

Following Ref. [30 let us make the following ansatz for 
the evanescent channel amplitude 



V'+(a;) ^}A„(j)nix) , 



(14) 



where the A^s are constants and the 0„'s are the bound 
state wave functions obtained by neglecting interchannel 
mixing in Eq. (I13p 



2mQ dx^ 



-£r 







(15) 



Notice that if Vm is also neglected in Eq. ([T^ the propa- 
gating channel corresponds to a free particle in ID. Then, 
in terms of the free-particle Green function, we may write 
the general solution of Eq. dE 



^- 



Akx 



ih k 



dx' 



[Vmt 



+ \x' 



(16) 



where k = ^/ImoE/h and [KkV'-hIk' denotes the action 
of the gradient-dependent potential on tp^ at point x' . 



Using now the ansatz (|14p in Eq. (|16|) , substituting in Eq. 
P^ and projecting on the set of bound states {(f>n,n = 
1, . . . , Nb} we obtain a matrix equation for the A„'s 



Nh 



Yl [('^"i ^ -^)'^«i«2 - A^nms] ^"2 = ^n, , (17) 

"2 = 1 



where 



mp 
ih k 



I dxidX2[Vm(l>nAxi[yr 



vaVn-i\x-i^ 



^ik\xi-X2\ 



and 



B„ 



dx [Vm4'ni]x e 



ikx 



(18) 



(19) 



Taking the limit x — > oo in Eq. 
tude of the transmitted wave 



16)) we find the ampli- 



t = 1 



mp 
ih k 



E^"^" 



(20) 



and the corresponding transmission T = |tp. The solid 
line in Fig. |4] lower panel displays numerical results ob- 
tained by solving the matrix equation p7)) , showing clear 
transmission minima when the energy is close to a bound 
state. A more explicit role of the bound states can be 
seen neglecting nondiagonal terms of the matrix Ai. In 
this case, the transmission amplitude reads 



t= 1 



mp Y^ 



\B„ 



ih k 



En-E -Mr 



(21) 



When E K, En the denominator in the right hand side of 
Eq. (j2ip reaches a minimum, thus yielding the mecha- 
nism by which the bound states produce deep minima in 
transmission. Notice also that A^„„ plays the role of a 
complex "self energy" that slightly distorts the position 
of the minima. Nevertheless, displacements of the dips 
from the bound state energies are hardly seen for weak 
Rashba couplings since Mnn ~ Q^o- I* *^^^ ^1^° ^^ shown 
that the relation Im(A^„„) — —mo\Bn\'^ /h^k is fulfilled 
and that this implies an exactly vanishing conductance 
at the dip energies. 

In the AP configuration no potential well explicitly ap- 
pears in the Hamiltonian, as shown is the lower panel of 
Fig. [TJ Nevertheless, the results of Fig. [5] prove that 
the xAP configuration also shows clear oscillations, with 
conductance dips in similar positions to the xP-polarized 
case. We can explain this quasibound states as a result 
of the combination of two effects: a) the refiection on the 
potential steps in vj^ and V- ; and b) the Rashba induced 
spin fiip. Indeed, adequately combined, the refiection 
and the spin flip may lead to a trapped state of the elec- 
tron. Mathematically, we could describe this mechanism 
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FIG. 5: (Color online) Conductance as a function of Rashba 
coupling intensity for polarization along a; in P and AP config- 
urations. Solid and open symbols are, respectively, the results 
with and without the Rashba mixing term a{x)pyax of Eq. 
([2|). Differently to the preceding figures, we take here the 
Fermi energy E as energy unit, with a corresponding length 
unit Lu = y'h'^/moE. The Rashba-coupling unit is then 
au ~ y'h'^E/mo. The remaining system parameters are: 
d = 20VELu, L = 8V5Lu, Ao = 20E. 
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(Color online) Same as Fig. [5] for polarization along 



by transforming the Hamiltonian with a local spin ro- 
tation 'D{x) = e~^'^''^^^\ where 0(a;) evolves from zero 
in the left contact to tt in the right one. In the trans- 
formed problem one component is then effectively bound 
by the two original potential steps. The transformation 
is rather cumbersome, however, due to the noncommu- 
tation of the kinetic term with (j){x), in addition to the 
also noncommuting Pauli matrices. 

In the y orientation of the contacts there is no coupling 
between -I- and — spins for q — 0, as immediately noticed 
from the term a{x)q{s\ax\s') of Eq. ([8]). This explains 
why there are no clear Fano oscillations for energies below 
threshold in Fig. [3] The minor features, can be attributed 
to the Rashba-induced coupling for finite g's. 



IV. DATTA-DAS TRANSISTOR 

The Datta-Das transistor— relies on the oscillatory 
character of the conductance as a function of the Rashba 
intensity ao- In this section we discuss the dependence on 
ao, considering first fully polarized contacts {E < 2Ao) 
and, subsequently, partial polarization at the end of 
the section. The fully polarized results agree overall 
with other theoretical analysis^^""— and, qualitatively, 
with the experiments of Koo et al.^ Quite surprisingly, 
however, the oscillatory character of the conductance is 
rapidly washed out if partial polarization is considered 
in our model by increasing the energy above the Zeeman 
threshold E > 2Ao. 

Figure [5] shows the results of our model for fully polar- 
ized leads with spin oriented along x. Upper and lower 
panels correspond, respectively, to xP and xAP configu- 
rations. In each case, solid symbols represent the results 
for the full Rashba Hamiltonian while open symbols cor- 
respond to the neglect of the mixing term. The con- 
ductance shows a damped sinusoidal behavior in both 
cases, with decreasing amplitude as ao rises. These re- 
sults agree with the expected Datta-Das behavior and, 
therefore, confirm the precession scenario in the contin- 
uum 2D case. Similar damped oscillations were obtained 
in Refs. [23 - l24l . Notice also that this damping is due to 
the Rashba mixing since it is absent in the open sym- 
bols. The oscillation period changes with the distance £ 
and successive minima approximately fulfill the spin pre- 
cession condition iao — nirfi /rriQ, with n = 1, 2, . . .. In 
AP configuration the conductance vanishes when ao = 
as a consequence of the spin mismatch between both con- 
tacts, known as spin-valve behavior. In the presence of 
Rashba coupling, however, the spin valve behavior is de- 
stroyed and we observe that the conductance rapidly in- 
creases with ao, at small couplings, and then oscillates is 
a similar way to the P case (Fig. [5] lower panel). 

For polarized leads along z the results are very sim- 
ilar to those already discussed for polarization along x 
and, thus, they will not be shown. On the contrary. Fig. 
ini corresponds to the configurations along y. Like be- 
fore, upper and lower panels are for yP and yAP con- 
figuration while solid and open symbols represent the 
cases with and without band mixing, respectively. In this 
case the conductance oscillations are almost absent, spe- 
cially in the yV arrangement (upper panel) , a result that 
agrees with the experiments* and with the precessing 
spin scenarioi^ The yP conductance decreases very slowly 
with ao and the effect of mixing is minimal, around a 10% 
decrease. In j/AP orientation we see how the Rashba mix- 
ing term again destroys the spin- valve effect at finite ao 's 
and we observe a small increment in conductance as the 
Rashba intensity increases. In this configuration there is 
a reminiscence of the oscillating behavior although much 
weaker as compared with the x ov z orientations. 

We turn now to the partially polarized cases, when the 
energy condition E > 2Ao allows both spins to propagate 
in the contacts. First notice that the polarization in a 
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FIG. 7: (Color online) Conductance as a function of ao for 
partial polarizations in a;P configuration. The different curves 
are for varying partial polarizations, from full (p = —1) to 
20% {p = -0.2). Parameters: d = 20\/5Lc7, L = 9>^/ZLu, 
units as in Fig. O 



given contact c — L^Ris given by 



Pc(^,Ae) = 



A. 



I A, 



(£;>2|A,|) 



{E < 2|A, 



(22) 



where, as already mentioned in Sec. II, we define 
(Ai,Aij;) to be (Ao,Ao) in the P configuration and 
(Ao,— Aq) in the AP configuration and Ao is assumed 
positive. The xP results for partial polarizations are pre- 
sented in Fig. [71 Notice that the oscillatory behavior 
is greatly quenched when the polarization is decreased, 
being heavily damped at \p\ = 0.5 and totally washed 
out at \p\ — 0.2. Thus, at low polarizations, our model 
predicts a monotonous decrease of the conductance with 
the intensity of the Rashba coupling that is not consis- 
tent with the operation of the Datta-Das device;^ This 
result shows the importance of having a high degree of 
polarization in the ferromagnetic contacts for obtaining 
a robust sinusoidal behavior. 

The above results are not substantially modified when 
using other system parameters, such as changing the 
energy or the distance between the leads d. The 
monotonous decrease of the conductance, without oscil- 
lations, at low polarizations is also seen in xAP, yP and 
yAP configurations. 



V. SPACE-DEPENDENT EFFECTIVE MASS 

In this section we investigate the relevance of having 
different effective masses in the semiconductor central re- 
gion and the polarized contacts. Till now the contacts 
were considered semiconductor materials with a Zeeman 
field in a given direction. A generalization towards fer- 
romagnetic materials in the contacts has to include the 
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FIG. 8: (Color online) Same as Fig. [2] assuming the effective 
mass in the contacts is given by rric = 15mo (c = L,R). 



different effective masses of ferromagnet and semiconduc- 
tor. As in Ref. |3l|, we then consider the effective mass 
in the contacs is the bare electron mass rue while in the 
central region it is given by the conduction band effec- 
tive mass of the semiconductor, 0.023TOe for InAs and 
0.063me for GaAs. Our aim is not a realistic modeling 
of ferromagnetic contacts, but to explore the qualitative 
effects of a position dependent mass on the preceding full 
semiconductor scenario. In particular, we shall still vary 
the energy from full to partial polarization, which is not 
a very realistic assumption for a ferromagnet. 

In the spirit of the model, we now use a generalized ki- 
netic term with a position dependent effective mass m{x) 
(see Appendix E| evolving from the semiconductor mass 
Too in the central region to 15mo in both contacts; 



n 



dx 2m{x) dx 2m{x) dy"^ 



(23) 



The big jump in effective mass at the interface is 
smoothed using Fermi functions as explained in Ap- 
pendix |X1 The presence of these effective-mass inter- 
faces is an additional source of conductance oscillations, 
as compared to the discussion of the preceding sections. 
Indeed, in this case even the ao = conductance with 
fully polarized contacts displays Ramsauer oscillations, 
as shown by the black symbols in the upper panel of Fig. 
m Based on the preceding section results, the addition of 
the spin orbit coupling is expected to introduce new oscil- 
lations of Fano type due to the coupling with quasibound 
states. Surprisingly, both types of oscillations interfere 
destructively, specially in the vicinity of the polarization 
threshold E — 2Ao, as shown by the gray (red) sym- 
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FIG. 9: (Color online) Same as Fig. [7] assuming the effective 
mass in the contacts is given by nic = 15mo (c = L,R). 



bols in Fig. |S] upper panel. Another conspicuous effect of 
the effective mass discontinuity is the big enhancement 
of conductance when the energy exceeds 2Ao. This is 
clearly noticed when comparing the upper panels of Figs 
[S]and[ll 

The lower panel of Fig. [S] shows the xAP conductance 
with position dependent effective mass. As a difference 
with Fig. [2] (lower panel), there are Ramsauer oscillations 
due to the mass jumps for E > 2Ao even for a vanishing 
ao- Below threshold we find Rashba- induced oscillations 
that look very similar to the xP ones in the upper panel. 
For y oriented contacts the results (not shown), as com- 
pared to those of Fig. [3l are also characterized by the 
appearance of clear Ramsauer oscillations below thresh- 
old in the yV configuration while in the j/AP orientation 
the variations are much smaller. 

A natural question to ask is whether the effective mass 
modification affects the conductance oscillations with ao 
discussed in Sec. IV. This is addressed in Fig. [9] for 
the fully polarized xP configuration. There are small 
changes, of course, but the overall behavior with damped 
oscillations is well preserved. Another result we should 
check is the disappearance of the oscillations at partial 
polarizations of the contacts (Fig. [7]). As proved by 
Fig. [51 this result is also robust with respect to effec- 
tive mass changes. Actually, the quenching of the oscil- 
lations at partial polarization is enhanced when the mass 
in the contacts is taken to be the bare mass: already for 
\p\ < 0.8 the conductance becomes monotonous, having 
only a slight decrease with ag. 



VI. CONCLUSIONS 

The linear conductance of a 2D stripe of Rashba inter- 
action with polarized contacts behaves in characteristic 
ways depending on the electron energy, the Rashba inten- 
sity as well as on the spin orientation and the degree of 
polarization of the contacts. We have explored all these 
degrees of freedom using a single model based on an ef- 



fective mass approach combined with phenomenological 
Rashba coupling and Zeeman fields. Our analysis was 
mostly numerical, in order to avoid additional simplify- 
ing assumptions. Nevertheless, semianalytical analysis 
in terms of coupled channels and quasibound states was 
also provided. 

We have identified two types of oscillations: Ram- 
sauer oscillations due to discontinuities in the effective 
potential and effective mass; and Fano oscillations due 
to the coupling with quasibound states. The latter ones 
are exclusively due to the Rashba coupling. For ener- 
gies below the full polarization threshold and without 
mass jumps at the interfaces we obtained pure Fano os- 
cillations in both parallel and antiparallel polarizations 
along X, while they are absent for polarization along y. 
With the addition of mass jumps at the interfaces the 
oscillations are quenched; a result that we attribute to 
a destructive interference between Fano and Ramsauer 
oscillations. Above the polarization threshold the con- 
ductance displays Ramsauer oscillations in most cases. 

Regarding the oscillations in conductance as a func- 
tion of Rashba coupling, our main result is the rapid 
quenching of the oscillations when the contacts are par- 
tially polarized. This oscillation quenching is even more 
pronounced when the effective mass increases in the con- 
tacts. The oscillations disappear for polarizations below 
20% and 80% for constant and position-dependent effec- 
tive mass, respectively. 
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Appendix A: Smooth transitions 

This appendix contains the precise mathematical forms 
of the smooth functions a{x), A{x) and m{x) defining the 
system Hamiltonian. We model the step-like character of 
these quantities using Fermi functions 



■^XQ^a ['^ } 



1 



1 + q{x-xo)/o 



(Al) 



where xq is the position of the step and a is giving the 
length around xq in which the transition takes place. Pre- 
cisely, it is 

a{x) = ao [J^t>/2,a{x) - T-i/2,a{x)] , (A2) 

A(a;) = ALJ-_rf/2,,(x) + Afl, [l - J"d/2,.(x)] ,(A3) 

+ mo [l'd/2,a{x) - J^-d/2.a{x)] ■ (A4) 

where the constants are: a^, Rashba intensity; A^, A^, 
Zeeman fields; m^, mg and mji, effective masses. Notice 



that, for the sake of simphcity, we assume a common 
value for a in Eqs. ([X2\\ . !(M\i and (K^ . 



Appendix B: Resolution method 

The numerical calculation of the linear conductance 
Eq. (fTO]) at a given energy E involves two steps. First, for 
a certain angle 9, or what is equivalent, a certain trans- 
verse momentum q, the coupled equations for -0^+ and 
ipq-, Eqs. ([HI), are solved to obtain the transmissions Ts's 
and T^ig. This is accomplished using the transmitting- 
boundary algorithm as in Ref. |l^. The calculation of 
Ref. Il9| was for quantum wires, with a confinement po- 
tential in the transverse direction, where the system of 
coupled equations was infinite and had to be truncated. 
The present case is, in this respect, simpler since only 



the two spin components of a given transverse momen- 
tum need to be considered. Nevertheless, the reader is 
addressed to Ref. [l3 for the technical details on how the 
differential equations with open boundary conditions are 
transformed into a linear system of equations. ^^ 

Once the transmissions at a fixed 9 are obtained, a 
second step of the calculation requires to integrate over 
the angle to calculate the linear conductance from Eq. 
(flUl) . This integral turns out to be somewhat delicate 
due to the presence of resonances as discussed in the 
quasi-analytical solution by ansatz of Sec. III.B. The 
6'-integration is then carried out using Gauss-Legendre 
quadratures with a certain set of abcissae and weights. 
To make sure that the integral is well converged we keep 
increasing the number of Gauss-Legendre points until a 
required accuracy is reached in a stable way. Typically, 
we require the error to be AG/Gq < 10"'^. 
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